arXiv: 1503.00061 vl [astro-ph.CO] 28 Feb 2015 


Mon. Not. R. Astron. Soc. 000,[!]{T7](2015) Printed 3 March 2015 (MN [AT E X style file v2.2) 


Separating weak lensing and intrinsic alignments using radio 
observations 

Lee Whittaker, Michael L. Brown, & Richard A. Battye 

Jodrell Bank Centre for Astrophysics, School of Physics and Astronomy, University of Manchester, Oxford Road, Manchester Ml 3 9PL 


3 March 2015 


ABSTRACT 

We discuss methods for performing weak lensing using radio observations to recover infor¬ 
mation about the intrinsic structural properties of the source galaxies. Radio surveys provide 
unique information that can benefit weak lensing studies, such as HI emission, which may be 
used to construct galaxy velocity maps, and polarized synchrotron radiation; both of which 
provide information about the unlensed galaxy and can be used to reduce galaxy shape noise 
and the contribution of intrinsic alignments. Using a proxy for the intrinsic position angle of 
an observed galaxy, we develop techniques for cleanly separating weak gravitational lensing 
signals from intrinsic alignment contamination in forthcoming radio surveys. Random errors 
on the intrinsic orientation estimates introduce biases into the shear and intrinsic alignment 
estimates. However, we show that these biases can be corrected for if the error distribution is 
accurately known. We demonstrate our methods using simulations, where we reconstruct the 
shear and intrinsic alignment auto and cross-power spectra in three overlapping redshift bins. 
We find that the intrinsic position angle information can be used to successfully reconstmct 
both the lensing and intrinsic alignment power spectra with negligible residual bias. 

Key words: methods: statistical - methods: analytical - cosmology: theory - weak gravita¬ 
tional lensing 


1 INTRODUCTION 

Gravitational lensing describes how light from a distant back¬ 
ground object is deflected due to the curvature of spacetime induced 
by the presence of a foreground matter and energy distribution. It 
is a rich phenomenon and can manifest itself in many ways. Two 
commonly studied regimes are “strong lensing” and “weak lens¬ 
ing”. Strong lensing is concerned with the distortion of a set of 
background galaxies or galaxy clusters due to a large foreground 
mass distribution close to the line of sight. This regime results in 
multiple images of galaxies (first detected with the observation of 
the Twin Quasar ( [Walsh et al,|1979| ) or, under the right conditions, 
Einstein rings. Weak lensing is concerned with the small but co¬ 
herent distortions of a set of background objects and, as such, is 
difficult to detect [Kaiser & Squires|1993[ l. 

The distortion of an image can be described by a linearized 
mapping that can be decomposed into two components - conver¬ 
gence and shear (see e.g. |Bartelmann & Schneider| ( j200l| ). On 
cosmological scales one is primarily interested in measuring the 
convergence power spectrum, which can be reconstmcted using es¬ 
timates of the shear. On these scales the shear is a result of the 
large-scale structure of the Universe and is known as cosmic shear. 
In 2000 cosmic shear was detected by four independent groups 


(Bacon et al. 2000 

Kaiser et al.|2000,|Van Waerbeke et al. 2000, 

Wittman et al. 2000 

). These detections laid the foundations for ob- 


servational weak lensing, which has been increasing in precision 
ever since and has now become a promising tool for probing cos¬ 
mology (jKilbinger et al.|20l3). 


The standard method for performing a cosmic shear measure¬ 
ment requires averaging over the observed ellipticities of a suffi¬ 
cient number of background galaxies and assuming that the average 
ellipticity is the consequence of cosmic shear. This method is built 
on the assumption that there is a zero intrinsic alignment (IA) of 
the source galaxies. However, non-zero IA signals were predicted 
as early as 2000 (Heavens et al. 2000 |Croft & Metzler|2000||Crit-| 


tenden et al. 2001 Catelan et al. 2001), with a subsequent detection 
being made in low-redshift surveys soon after this [Brown et al.| 
|2002[ >. The presence of IA has the effect of producing a false shear 
signal and hence a bias in the standard method. Although it has not 
had a significant impact on the present generation of surveys (e.g. 
Heymans et alTj ( [2012^ ), it is likely to be important in forthcoming 
surveys such as those performed by, for example, DES [The Dark 


Energy Survey Collaboration|2005), Eucli d fCimatti & Scaramella 
2012 1 and the LSST |LSST Dark Energy Science Collaboration 


2012i. 


Most weak lensing surveys so far have been performed in the 
optical waveband, but the SKA promises the possibility of perform¬ 
ing surveys in the radio waveband (Brown et al. |2015[ l. Such a 
survey offers some unique advantages. Firstly, based on an idea 
originally proposed by [Blain] ( |2002[ >, |Morales| [2006} introduced 
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a method for performing weak lensing using resolved galaxy ve¬ 
locity maps. Radio HI emission is the most promising part of the 
electromagnetic spectrum to construct the velocity maps, due to 
the brightness of the emission lines and the well understood lumi¬ 
nosity characteristics. Gravitational lensing leads to a velocity map 
which is inconsistent with the observed galaxy image and |Morales| 
< |2006) > showed that this effect can be used to recover estimates of 
the underlying shear signal. In principle, this method removes the 
contribution of both galaxy shape noise and the effects of intrinsic 
alignments from weak lensing surveys. It does, however, require 
velocity maps from well resolved galaxy images, which reduces 
the number density of available galaxies in the survey. However, 
using a toy model [Morales| i {20()6} showed that this method may be 
competitive in future radio surveys, such as with the SKA. 

Secondly, [Brown & Battyej < f201 lb| ) (hereafter BB11) sug¬ 
gested a new technique in radio weak lensing which would use the 
polarization information contained in the radio emission of a source 
galaxy as a tracer for the intrinsic position angle (IPA) of the galaxy. 
It had previously been shown ( [Dyer & Shaver|1992fr that the net po¬ 
larization position angle is unaffected by a gravitational lens. For 
future deep radio surveys, the population of observed galaxies is 
expected to be dominated by star forming galaxies. The dominant 
source of radio emission from such a galaxy is expected to be syn¬ 
chrotron radiation emitted as electrons are accelerated by the large 
scale magnetic fields within that galaxy. This emission gives rise to 
a net polarization position angle (PPA) which, on average, is anti¬ 
aligned with the plane of the galaxy fStil et al.|2009j l, providing in¬ 
formation about the galaxy’s intrinsic orientation. It was shown that 
such information can be used to construct a shear estimator which 
greatly reduces the biases resulting from intrinsic alignments com¬ 
pared to the standard method, and also reduces the errors on the 
shear estimates. It was shown i jBrown & Battye||201 la| l that this 
new method can be used to create foreground mass reconstructions 
with accuracies comparable with the standard method, subject to 
specific assumptions on the size of the errors on the estimates of 
the intrinsic orientations of the galaxies, and the fraction of galax¬ 
ies with reliable polarization information. In principle, the method 
can be applied to estimates of the IPA from any source, and a sim¬ 
ilar analysis could also be applied to the technique described by 
[Moralesjf2006| . 

The method displays great promise. However, there is a small 
residual bias in the estimator when there is both a non-zero error 
in the IPA estimates and a non zero IA signal. In this paper we 
develop improved estimators which remove this bias. In Section 
[2] we present an overview of the method proposed by BB11. We 
discuss the noise properties of the method and the residual bias 
which is introduced in the presence of both an IA signal and an 
error on the intrinsic position angle estimates. We address this bias 
by introducing a correction term, which depends on the form of the 
error distribution. 

In Section [3] w e extend the angle-only estimator, introduced 
by |Whi ttaker et al. (2014), to include IPA information and we also 
introduce a hybrid method that combines an angle-only estimate of 
the intrinsic alignment with full ellipticity information. In Section 
[4]we test the methods using simulations by reconstructing the shear 
and IA auto and cross-power spectra in three overlapping redshift 
bins. We conclude in Section[5] 

Throughout the paper we assume that we have reliable IPA 
information for every galaxy for which we have reliable ellipticity 
or observed position angle measurements. For a real radio survey 


this would not be the case. However, the purpose of this paper is to 
demonstrate the potential of the methods presented, provided that 
we have a sufficient number of galaxies to recover a reliable shear 
estimate. A detailed discussion of the fraction of galaxies expected 
to have reliable polarization information can be found in BB11. 


2 METHODS AND TECHNIQUES 

Working well within the weak lensing regime we can express the 
observed ellipticity of a galaxy, e° , as the sum of the intrinsic 
ellipticity, e lnt , the shear, 7 , and a measurement error, 5e , such that 

e =e +7 + <5e. (1) 

Considering a small region of the sky, or a cell, such that the shear 
can be considered constant within that cell, and assuming that the 
mean intrinsic ellipticity of the galaxies in that particular cell is 
zero, we can recover an unbiased estimate of the shear by averaging 
over the observed galaxy ellipticities: 

1 N 

t=s (2) 

i=l 

If the unlensed galaxies within a particular region of sky are 
randomly orientated, then the average intrinsic ellipticity is zero. 


However, there is theoretical motivation (Catelan et al. 200T][CriF 


tenden et al. 2001, ling 2002, Mackey et al. 2002, Hirata & Sel- 

jak 2004 

and observational evidence (Brown et al. 2002, Hey- 

mans et a 

. 2004, Mandelbaum et al. 2006,2009, Hirata et al.|2007; 


Brainerd et al.|2009) indicating that during galaxy formation, cor¬ 
relations between the intrinsic position angles of galaxies, a lnt , 
may arise if those galaxies share an evolutionary history. This is 
the source of the IA signal. If we now assume that the mean intrin¬ 
sic ellipticity (or, equivalently, that the IA signal) is non-zero, then 
the standard estimator of equation 0 is biased, such that 

<7-7> = (e mt ) • 0) 

If we assume that the intrinsic ellipticity of a single galaxy can be 
expressed as the sum of the intrinsic alignment signal, 7 IA , and a 
randomly orientated ellipticity, e ran , then 

int IA , ran / a \ 

e = 7 + e , (4) 

and the bias in the standard estimator becomes 

< 7 - 7 )= 7 IA - 0 ) 

Hence, an estimate of the shear recovered using the standard 
method yields a result which is biased by the IA signal. 


2.1 The Brown & Battye (BB) Estimator 

In order to mitigate the effects of the bias introduced by IA, BB 11 
proposed using polarization information from radio surveys to re¬ 
cover an estimate of the IPA, a mt . It was found that, for the ideal 
case where there is a zero error on the IPA measurement, and where 
the PPA is an exact tracer of a mt , the shear can be recovered ex¬ 
actly, using only two source galaxies. 

Expressing the intrinsic ellipticity in polar coordinates, the 
components of the observed ellipticity can be written as 

ef s = ]e int | cos (2a int ) + 7l + for, 

= je mt | sin ^2a lnt ^ + 72 + < 5 e2. (6) 
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Figure 1. The residual bias in the BB estimator from 10 4 realizations, with 10 4 galaxies in each realization. For each realization 71 , 
72 and 7 ^ a are selected randomly with a range [—0.1, 0.1]. The left panel shows the bias in the BB estimator as a function of cr a int, 
with 7 ^ a = 0.05. The right panel shows the bias as a function of 7 j A with (Ji nt = 15°. In both cases the red curve is the linear 
approximation of the bias, given in equation CD- 


If we define the pseudo-vector 


rii 


sin( 2 aj nt ) \ 
-cos( 2 a‘ nt ) )’ 


(V) 


where d mt is an estimate of the intrinsic position angle provided 
by a measurement of the IPA, a new estimator for the shear can be 
derived, such that 


7 = A 1 b, 


( 8 ) 


where A is a 2 x 2 matrix and 6 is a two-component vector 

N 

A = '^2w i n i nJ', (9) 

i= 1 


b = ^2wi (e° bs • hi ) hi, (10) 

2=1 

and wt is a normalized arbitrary weight assigned to each galaxy. 

In the presence of a non-zero IA signal and a non-zero error 
on the estimate of a , it is found that the estimator given in equa¬ 
tion § is biased, although this bias is suppressed significantly with 
respect to that of the standard estimator. If we assume that the com¬ 
ponents of the intrinsic ellipticity are isotropically distributed about 
the IA vector, 7 IA , and use a uniform weighting, such that w l — 1 , 
it is possible to gain some insight into the nature of this bias. If 
we make the further assumptions that IV 3> 1 and the IA signal is 
much smaller than the spread in intrinsic ellipticities, <r e , that is if 
| -y | IA -C <r e , we can approximate A to leading order in "y IA as 

N 

A « —I, (11) 

and hence the estimator can then be approximated as 

7 W Jf h ’ (12) 


The noise properties inherent in using measurements of the PPA 
as a tracer of a mt are discussed in BB11. For this discussion we 
assume that the measurement error, <5a mt , is independent of the 
true IPA and distributed symmetrically about zero. If we then make 
the substitution & lnt = a lnt + <5a lnt , we can write the expectation 
value of the trigonometric functions of the IPA as jWhittaker et al.| 
[20l4| 

(cos ( 2 a lnt j ) = (cos ( 2 a lnt ) ) ftp, 

(sin ( 2 a lnt j ) = (sin ( 2 a mt ) ) dT, (13) 

where 

/3n nt = (cos (n<5a mt ) ) , (14) 

which is the mean cosine of the distribution of <5a lnt , and where n is 
an integer. For a Gaussian measurement error, this can be simplified 
to 

= exp , (15) 

where o^int is the standard deviation of the measurement error and 
is expressed in radians. Taking the limit N —> 00 and using the 
result of equation ( | 1 3| >, it can be shown that equation ( | 12 | > may be 
expanded to first order in the shear and IA, such that the bias in the 
estimator is 

(7 - 7> « (l - $! nt ) 7 IA - (16) 

For a a int <C 1, one finds that (7 — 7 ) ~ 2 cr^ int 7 IA and, there¬ 
fore, we see that the bias is suppressed by a factor of 2a^ int relative 
to the standard estimator. The bias in the BB estimator is illustrated 
in Figure |T] where we assume a Gaussian measurement error on 
the estimate of a lnt and a Rayleigh distribution for the intrinsic el¬ 
lipticity distribution (which we define as the distribution of |e ran |). 
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such that 


In Appendix A it is shown that the matrix M* is given by 


nwn) = 


1 — exp I — 


2a? 


exp 


2 O'? 


(17) 

where |e™xl is the maximum allowed value of the modulus of 
e ran ; for all of the simulations in this paper we have assumed a 
Rayleigh distribution for 


with values of 


= 1 and 


a e = 0.3/\/2- From Figure [T] we see that the estimator success¬ 
fully reduces the bias introduced by the IA. There is, however, a 
residual bias introduced when both the measurement error on a mt 
and the IA signal are non-zero. 

In the limit | 'y IA | <C <r e the standard error for the shear esti¬ 
mator can be written as 


2a, 


2 (1 - Pt) 

N 


(18) 


where a is the measurement error on the components of e obs . As¬ 
suming this error to be zero, and assuming a a i D t -C 1 , the error 
can be approximated as 

_ ~ 4<Ta int C e int nm 

(19) 

and hence we see that 07 is suppressed by a factor of 4cr a i n t relative 
to the standard estimator, in agreement with the findings of BB11. 

Given an estimate of the shear, and assuming that the effects 
of the intrinsic ellipticity can be modeled using equation an 
estimate of the IA signal can be recovered trivially, such that 


- IA 

7 = 


1 N 


7- 


( 20 ) 


( 21 ) 


The bias in the estimate of the IA signal arises from the bias in the 
shear estimator and hence to first order, has the same magnitude as 
the bias given in equation but with the opposite sign. 

To first order in the shear and intrinsic alignment, the error on 
the IA estimator is due to the random shape noise, 

a e 

a~ ia = >—=. 

1 Vn 

From equation ( |2 1 1 ) we see that the error on the IA estimator is in¬ 
dependent of the error on a lnt and therefore, there is no suppression 
of this error by < 7 a int. 


2.2 The Corrected BB (CBB) estimator 

It is possible to construct an unbiased shear estimator in the limit 
N —> oo by following the approach outlined in BB11. This cor¬ 
rected form of the BB estimator (hereafter the CBB estimator) can 
be written as 


7 = D“ 1 h, 


where D is a 2 x 2 matrix 


D = ^M„ 


and where h is a 2 -component vector 

N 

h = J2 M l£ f s . 


( 22 ) 


(23) 


(24) 


M, = 


p?' - cos (4af‘) -sin (4A;' , 


v int ) 


— sin (4 a'i 


/3i nt + cos (4aj nt ) 


( 25 ) 


where the term ft™ 1 is defined in equation © and corrects for the 
bias on the trigonometric functions introduced by a measurement 
error on a lnt . Once one has an estimate of 7 , an estimate of the IA 
can be recovered using equation ( | 20 | ). 

We have tested the performance of the CBB estimator using 
simulations composed of 100 galaxies and assuming an input shear 
signal of 71 = —0.03 and 72 = 0.04, and an input IA signal of 
7 i A = 0.01 and 72 A = —0.02. We recovered shear and IA esti¬ 
mates from 10 4 realizations using the original form of the BB es¬ 
timator (equations © and ( |20| )) and the CBB estimator (equations 
(|20|l and 1 22 1 ). We assumed a zero error on measurements of e obs 


and a Gaussian measurement error with r.m.s. 10° on a ln . The re¬ 
sults of this test are shown in Figure [2] Table [T] presents the mean 
recovered shear and IA estimates and the standard deviation of the 
estimates. 

Figure [3] shows the residual bias in the CBB estimator. The 
bias correction term, /34 nt , corrects for the bias introduced to the 
mean trigonometric functions in the BB estimator when there is 
an error on the estimates of a mt . However, for a finite number of 
source galaxies, there will also be noise in the estimates of the mean 
trigonometric functions (which enters into the CBB estimator via 
the inverse of matrix D). This noise propagates nonlinearly into 
estimates of the shear, resulting in a residual bias which is not cor¬ 
rected for. In the tests we have conducted, we find that this residual 
bias is always much smaller than the dispersion in the shear es¬ 
timates and contributes to a negligible residual bias in the power 
spectra reconstructions discussed in Section[4] 

In order to estimate the dispersion in the shear and IA esti¬ 
mates, we can write an approximate form of the CBB estimator. To 
leading order in 7 and 7 14 , the CBB estimator can be written as 


7i = 


NP^ 


E 


.oint obs,(i) obs,(i) 
P 4 e l ~ e l ( 


(4a' nt ) 


s 5 ) • (a ~ int \ 

sin ( 4 OLi J 


1 

^ =N0? 


, obs,(i) . 

+ e 2 cos I 


^nt e ob s ,(i) _ £ obs,(») gin 

(4ai nt )j, 

1 N 

- IA -L \ ^ T obs,(i) ( a ~int\ . obs,(i) . ( A ~int\l 

71 C0S J +e 2 sm ( 4 «i JJ . 

i= 1 

1 N 

~IA ^ \ A [ obs,(i) . ( . ^int\ obs,(i) / . ~int\l 

72 =^L[ £ 1 sm [4tti )-e 2 ^ cos (4a, Jj. 

^ i=1 

(26) 

The error on the CBB estimator can then be approximated as 

(l-PT t2 ) +a 2 (l + /r 2 )' 

Np\ nt2 


( j^ia =a*,iA = 
>i >2 


2 , 2 

(7 e (T 


(27) 
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Figure 2. The recovered shear and IA estimates from 10 4 realizations, with each realization consisting of 100 galaxies, and assuming 
a measurement error on the IPA of tr i n t = 10°. The black curves show the distributions of recovered shear and IA estimates when 
using the CBB estimator (equations (20) and (22) ) and the vertical black lines show the mean recovered shear estimates when using this 
estimator. The red curves show the distributions of recovered shear estimates when using the original BB estimator (equations {8) and 
(20) ) and the vertical red line shows the mean recovered shear estimates when using this estimator. The green dashed lines, which lie on 
top of the black lines, show the input shear signal. The success of the correction to the original BB estimator is clearly visible in these 
plots. There is, however, a modest increase (~20%) in the dispersion of the shear estimates and a ~30% increase in the dispersion of 
the IA estimates when using the corrected form of the estimator with this set of input values; this is quantified in Table]!] 


Estimator (x 10 2 ) 

°71 

^12 

<71 > 

(72} 

a i\ K 


<7i A > 

<72 A > 

Original BB 

1.40 

1.40 

-2.79 ±0.01 

3.58 ±0.01 

2.13 

2.12 

0.71 ± 0.02 

-1.55 ±0.02 

Corrected BB 

1.71 

1.70 

-3.00 ± 0.02 

4.02 ± 0.02 

2.75 

2.73 

0.92 ±0.03 

-1.99 ±0.03 


Table 1. The mean and standard deviation of the shear and IA estimates recovered from 10 4 simulations. Values are quoted for both 
the original BB estimator (equations {8) and (20) ) and the CBB estimator (equations (20) ) and m The input shear and IA values are 
71 = -0.03,72 = 0.04, 7 } a = 0.01 and 7^" = -0.02. 
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Figure 3. Same as for Figure[T]but for the CBB shear estimator. We see a residual bias which is a result of the finite number of source 
galaxies, however, this residual bias is much smaller than the residual bias in the original BB estimator, shown in Figure]!] 



Using equation {27} with the input values used to produce Fig¬ 
ure [2] we recovered approximate values for the dispersion in the 
estimators of a *. = a*, 0 = 1.68 x 10 -2 and cr- ia = a-iA = 

II iz 7 1 7 2 

2.71 x 10 -2 , which are in good agreement with the measured 
values quoted in Table [T] For completeness, we also used equa¬ 
tions l |18[ > and ( |2 1 1 > to recover approximate values for the disper¬ 
sion in the original BB estimator. These values were a* ll = cr^ 2 = 
1.40 x 10~ 2 and <t-ia = <r- ia = 2.12 x 1CP 2 , which are also in 

“i 72 i—i 

agreement with the values in quoted TablellJ 


2.3 Required Galaxy Numbers for the CBB Estimator 


The CBB estimator becomes unstable when there is a low number 
of background galaxies available in a particular cell. To gain some 
insight into the source of this issue, we can examine the behaviour 
of the determinant of matrix D when a low number of galaxies is 
considered. The determinant of matrix D is 


det (D) =/?r 2 


1 N 
sE' 


(4<S,“) 


1 

N 


sin ^4di lnt ^ 


(28) 


As the measurement error on a mt is increased, the bias correction, 
Pi, decreases. For a finite number of background galaxies, chance 
alignments of the random components of the intrinsic galaxy orien¬ 
tations can force this determinant to approach zero, with the effect 
being more likely when the number of galaxies in a cell is low. This, 
in turn, can produce substantial outliers in the estimated shear val¬ 
ues, as the modulus of a particular shear estimate is scaled with 
the reciprocal of the determinant. It is possible to place constraints 
on the number of background galaxies required for a reliable shear 
estimate by assuming that there are enough galaxies in the sam¬ 
ple, such that the central limit theorem can be applied to the dis¬ 
tributions of the means of the trigonometric functions in equation 
{28}. We can use this assumption to examine the probability that 


the sum of the square of the mean trigonometric functions in equa- 
tion 1 28 1 will lie within a given range of . The determinant 
is independent of the shear and is only dependent on the IA sig¬ 
nal at 4 th order, so we can safely assume the IA signal to be zero. 
With these assumptions in place, we choose to constrain the num¬ 
ber of galaxies, such that we can exclude values of the reciprocal 
>2//? 4 nt at a confidence level of 99.99994%, which is equivalent 
to a confidence level of 5a for the Gaussian distribution. The choice 
of 5cr is selected to mitigate the issue of outliers when considering 
the simulations in Section [4] where we reconstruct the shear and 
IA auto and cross-power spectra using ~10 G cells per redshift bin 
and therefore expect typically one cell per reconstruction to have 
a reciprocal value > 2/(5^ . The constraint on the values of the 
reciprocal >2//34 n is somewhat arbitrary but serves to provide 
an upper limit on the dispersion of the shear estimates. 


This choice of constraint parameters results in Figure]?] where 
we plot the number of galaxies required in the sample as a function 
of <T a int. As an example, let us assume a measurement error on 
a mt of 10°. Then, from Figure]?] we find that we need ~46 galax¬ 
ies in each cell, such that values of the reciprocal of the determinant 
>2//? 4 nt are ruled out at a confidence level of 5a. For Figure^we 
considered 100 galaxies per realization and hence outliers were not 
an issue for these tests. The number density of background galax¬ 
ies will be fixed for any specific set of observations. However, for a 
fixed number density of galaxies, the size of the cells may be cho¬ 
sen so that the number of source galaxies within each cell is greater 
or equal to the number of galaxies required to recover a reliable 
shear estimate. For a low number density of background galaxies, 
this will of course result in a large cell size and hence the loss of 
small scale information. 
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Figure 4. The number of galaxies in the sample, as a function of the er¬ 
ror on a int , such that the reciprocal of the determinant < 2//3J nt with a 
confidence level of 5cr. 


3 ALTERNATIVE APPROACHES 


We can also recover an estimate of F\ (| ~/ lA |) > 



( N \ 

-A 

> 

II 

^1- 

J2 cos ( 2d i nt ) 

\i=l J 


+ ( 5Z Sin ( 2 ®i nt ) 


(32) 

which can be inverted to provide an estimate of 17 IA |. Using equa¬ 
tions (|31| i and ( |32[ > we can, therefore, recover a complete estimate 
of 7 IA using only measurements of the IPA. 

The Fi (17 IA |) function describes the mean cosine of the an¬ 
gle between the vectors e ran and 7 IA as a function of |7 IA |. If we 
assume that the effects of IA can be modeled using equation 0- 
then the F] (17 IA |) function is found to be 


/I ta l\ 1 /“l e max| r 2 

Fl \r\)nJ J dQ d l eran l/(|e ran l) 




(\ IA 1 1 



X 

9 1 (J7 

| ! I 6 

J , 

where the function g 1 

(|V A |, 

|e ran | 

, a ran ) is given ; 


(\ ia| 



) ^ 


91 v 7 

.l e 


1 s/t'i+t'2 

with 






t 

€l = 

ia| . 
7 + 

1 ran i 

|e 1 

cos (2a ran ), 


(33) 


(34) 


£2 =|e 


| • /r» r 

| sin (2a 


(35) 


3.1 Full Angle-Only Estimator (FAO) 

In this section we extend the angle-only shear estimator, introduced 
by |Whittaker et al.| ( |2014[ >, to include measurements of the IPA. 
Assuming a prior knowledge of the intrinsic ellipticity distribution, 
|Whittaker et al.||2014| l showed that it is possible to recover an esti¬ 
mate of the shear using only measurements of galaxy position an¬ 
gles. Using measurements of the IPA, as opposed to the observed 
position angles, this method can be used to recover a direct estimate 
of the IA signal. 

We begin by writing the IA in polar form, such that 


IA 

7i 


IA 

72 


IA 

7 cos 


IA, 

7 | sm 


K) ■ 

in (2 q ia ^ . 


(29) 


It can then be shown that the means of the cosines and sines of the 
intrinsic position angles can be written as 

cos (2a lnt j^ =Fi (|7 IA |) cos ( 2a IA ^ ’ 
in (2a int ^ =F 1 (|-y IA |) sin (2a IA ) , (30) 

where the form of the function Fi (17 IA |) is dependent on 
/(|e ran |) and the assumed model which describes how the IA 
transforms e ran ->■ e int . Note that Fi (|7 IA |) is a function of the 
modulus of the IA only. Assuming a zero measurement error on the 
intrinsic position angle estimate (i.e. d lnt = a lnt ) and a sample of 


such that 


d IA = - tan 1 
2 


Ef= 


sin (2d‘ nt ) 


E i= i cos (2d‘ nt ) 


If we now allow for a measurement error on a mt , such that 

+ <5a in \ (36) 


where <5a lnt is independent of the true o' 111 , |Whittaker et al. ' 20141 
show that estimates of a IA remain unbiased. However, estimates of 
7 IA |, obtained by inverting the Fi (7 IA ) function, become biased. 
The bias can be corrected for by dividing the Fi (17 IA |) function 
by the correction term p ?'*'> which follows the definition given in 
equation (| 1 4|) : 




(M) 


N/3 “ 


1 


/E cos (2d( nt ) J -f I sin (2d 


(37) 

Assuming that we are working well within the weak lensing 
regime, such that e obs can be described using equation (JTJ, and ig¬ 
noring measurement errors, we can express the observed ellipticity 
in terms of the shear and IA, such that 

obs i IA , ran /oo\ 

e =7 + 7 + e . (38) 

Expressing the observed ellipticity in polar coordinates: 

obs 

obs I obs I • (n obs \ 

e 2 = e sm 12a 1 , (39) 

we can follow the approach described above to recover estimates 
of the vector 7 + 7 IA from the observed galaxy orientations. If we 
assume a measurement error on a obs which is independent of the 
true value 


obs 

6 

COS 

(2a obs ) 

obs 


^2Q obs ) 

6 

sm ( 



z' obs obs , c- obs 

a = a. + da , 

(40) 


we can define the terms /3£ bs , such that 


(31) 

/?° bs = /cos | 'nSa° bs ^\ \ . 

(41) 
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Upon expressing the vector 7 + 7 IA as 


to obtain the parameter, u, analytically. This is found to be 


, IA I tot I /o tot\ 

7i+7i = 7 | cos (2a j , 


f 7T \ 

, IA I tot I • /o tot\ 

72 + 72 = 7 sin (2a j, 

(42) 



( 47 ) 


we can recover an estimate of a tot , such that 


/'tot 1 , —1 

a — tan 
2 


Ef =1 sin(2a° b3 ) 
Eticos (2d° bs ) 


(43) 


and an estimate of 7 tot which satisfies the equation 


* (77) = 


NP § 


Ml§' 


(2a° bs )J + (j>in (2d? bs ) 

(44) 

and which provides us with an estimate of the vector 7*°* using 
measurements of the observed position angle, d° a , only. 

The vector e ran , given in equation ( |38^ , is identical to that 
given in equation 0, therefore the form of the F\ (|7 tot |) 
function, given in equation (|44|), is identical to the form of 


the F\ (| 7 ia |) function in equation 133 1, with the substitution 
7 IA | ~+ Tlle term hi** corrects for the bias introduced 


mates of 7* A and 7 tot , an estimate of the shear can be recovered 
trivially, such that 

v IA 


7 : 


^ tot 

7 -7 


(45) 


To summarize, the full angle-only estimator (hereafter the 
FAO estimator) first requires an estimate of the intrinsic ellipticity 
distribution, / (|e ran |). We can use this information with measure¬ 
ments of the IPA only to recover an estimate of 7 IA via equations 
< |3 1 ^ and l |37| >. An estimate of the vector 7 + 7 IA can also be ob¬ 
tained using the same method via equations ( |43| ) and ( |44[ >. Finally, 
we use equation to recover an estimate of 7. 

Assuming the same input values as used in Figure[2] we recov¬ 
ered shear and IA estimates from 10 4 realizations using the FAO 
estimator. The error on a mt was assumed to be 10° and the error 
on q° s was assumed to be zero to allow for a direct comparison 
of the performance of this estimator with the CBB estimator, where 
we assumed zero errors on the ellipticity measurements (e obs ). The 
results of this test are shown in Figure[5] Note that the reduction in 
the dispersion of the shear and IA estimates is a result of the fact 
that we have assumed a perfect knowledge of the intrinsic ellipticity 
distribution. Errors on the prior knowledge of the intrinsic elliptic¬ 
ity distribution introduce multiplicative biases to the estimates of 
the shear and IA. This issue is addressed in |Whittaker et al.| ( j2014| , 
where constraints are placed on the size of the errors on the mea¬ 
surements of the ellipticities and the size of the sample used to 
estimate the intrinsic ellipticity distribution, such that this multi¬ 
plicative bias is below a desired threshold value. Table[2]shows the 
mean recovered shear and IA estimates and the standard deviation 
of the estimates. 

A linear form of the estimator can be obtained by following 
the approach outlined in Whittaker et al. (2014). Assuming that the 
Fi (I 7 I) function can be approximated as 

-Pi (|t|) « «|7l, (46) 

for a general intrinsic ellipticity distribution we can find the coef¬ 
ficient, u, numerically. However, assuming a Rayleigh distribution 
for | e ran | and assuming that cr e is small enough for us to safely al¬ 
low the limit in the integral, | eJJJlx | > to tend to infinity, it is possible 


The linear approximation of the FAO estimator can then be written 
as 


„ 

1 

N 

'cos (2d° bs ) 

cos (2d( nt ) 

7l 

uN 

2^ 

i=1 

PT S 

PT 


1 

N 

'sin (2d° bs ) 

sin (2d( nt )' 

72 

uN 

i=1 

Tq 

too 

of 

PT 

-IA 

7i 

1 

“uTV 

cos (2a( nt ) 

/ ^ flint ’ 

i=l 


-IA 

72 

1 

“uTV 

sin (2d‘ nt ) 

/ ^ flint 

i=l 



From here it is possible to recover an approximation for the disper¬ 
sion in the estimator, given by 


71 uVN [2P° hs2 

1 


/"T2 


+ 


+ 


~uVN 1 2/3§ bs2 2$, nt2 


1 

c* IA =ct-ia = —= 

7l 72 uVN 


2/)- a 

1 

;ir 

1 


APT 


1 


2 

^cos 

(2a obs ( 

cos 

(2a int )) 

2 






1 

2 

(sin 

( 2 a ob + 

sin 




(49) 


The dispersion in the shear estimates depends on the correlations 
between the cosines and sines of the true observed and intrinsic 
position angles. These, in turn, depend upon the intrinsic elliptic¬ 
ity distribution, the IA and the shear. In the absence of a shear 
signal, these correlation terms will, to first order in the IA sig¬ 
nal, equal 1/2. If we also neglect measurement errors, such that 
PT S = PT = 1, then the dispersion in the shear estimates be¬ 
comes zero. However, the presence of a non-zero shear signal re¬ 
duces the correlation between the trigonometric functions and an 
error is introduced to the estimates. Hence, in the absence of mea¬ 
surement errors on the position angle measurements, the leading 
order term in the dispersion is dependent on the true shear. Mea¬ 
surement errors on the position angles also increase the dispersion 
in the estimates, as expected. Using the input values assumed in 
Figure [5] with equation {49}, we recovered approximations for the 
errors on the shear and IA estimates by calculating the correlation 
terms numerically. The errors were found to be cr 71 = 1.06 x 1CP 2 , 
ox, = 1.13 x 10~ 2 and cr ta = cr ta = 2.54 x 1CU 2 , which are 

7l 72 1—1 

in good agreement with the values quoted in Table HI 

Assuming that the shear signal is zero, we can also recover 
estimates of the dispersion in the shear estimates, which are cr 71 = 
0.87 x 10~ 2 and cr 72 = 0.85 x 10 -2 . These values are approx¬ 
imately 20% lower than the values quoted above where shear was 
included. Hence, we can conclude that the dispersion in the shear 
estimates depends strongly on the input shear signal, even if mea¬ 
surement errors on the position angles are included. This is an issue 
when trying to remove noise bias in power spectra estimates, and 
is discussed in more detail in Section [4] The dispersion in the IA 
estimates are independent of the true IA signal to first order. 

Figure[6]shows the residual bias in the FAO estimator. As with 
the CBB estimator, it is expected that there will be some residual 
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Figure 5. The recovered shear and IA estimates from 10 4 realizations, with each realizations consisting of 100 galaxies, and assuming 
a measurement error on a int of cr a i n t = 10°. The black curves show the distributions of recovered shear and IA estimates when using 
the FAO estimator. The vertical black lines show the mean recovered estimates using this method. The red curves show the distributed 
shear estimates when using the hybrid method, with the IA estimates identical for both methods. The vertical red lines show the mean 
recovered estimates using this method. The dashed green lines show the input signal. Here we see that both methods have successfully 
recovered shear and IA estimates with negligible bias. The dispersion in the FAO shear estimates is ~40% lower than those recovered 
using the CBB estimator and the dispersion in the IA estimates is ~6% lower for this set of input values. The dispersion in the hybrid 
shear estimates is ~15% lower than those recovered using the CBB estimator. The results are quantified in Table[2] 


Estimator (x 10 2 ) 

ct 71 

a i'2 

(71 > 

( 72 } 

a i\ K 


<7i A > 

<7j A > 

Angle-only 

1.02 

1.07 

-3.02 ± 0.01 

4.02 ±0.01 

2.56 

2.57 

0.99 ± 0.03 

-1.99 ±0.03 

Hybrid 

1.42 

1.42 

-3.00 ± 0.01 

4.01 ± 0.01 

2.56 

2.57 

0.99 ±0.03 

-1.99 ±0.03 


Table 2. The mean and standard deviation of the shear and IA estimates recovered from 10 4 simulations. Values are quoted for both the 

S 'e-only estimator (equations {31} and {37}, and equations {43} - {45}) and the hybrid estimator (equations (H) and HZ)> and equation 
). The input shear and IA values are 71 = —0.03, 72 = 0.04, 7^ = 0.01 and 7^ A = —0.02. 
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Figure 6. Same as for Figure|T|hut for the FAO shear estimator. From this plot see that any residual bias in the estimator can be considered 
negligible. 


bias from the nonlinear propagation of noise on the mean trigono¬ 
metric functions into estimates of the shear. However, in all of the 
tests conducted this residual bias is found to be negligible. 


3.2 Hybrid Method 

In this subsection we introduce a hybrid method combining the 
standard method, which averages over galaxy ellipticity measure¬ 
ments, with the angle-only IA estimator. 

Using a knowledge of the intrinsic ellipticity distribution, 
/ (I 6 ™I)’ with measurements of the IPA only, we first recover an 
estimate of the IA signal via equations ( |3 1 1 > and ( |37| ). We can then 
combine this estimate of the IA signal with an estimate of the vec¬ 
tor 7 + 7 IA , provided by the mean of the observed ellipticities, to 
recover an estimate of the shear: 

7=(^E c ? bS )-^ A (50) 

Using the same set of realizations as used to test the FAO 
method in Figure [5] (black curves), we recovered shear estimates 
from 10 4 realizations using the hybrid shear estimator (equation 
(|50[). The error on a mt (which can be estimated using a measure¬ 
ment of the PPA) was assumed to be 10° and the error on e obs was 
assumed to be zero. The results of this test are shown in Figure [5] 
as the red curves. It should be noted that, since the same realiza¬ 
tions have been used to test the hybrid and FAO methods, the IA 
estimates are identical. As for the FAO estimator, discussed in Sub¬ 
section ED the reduction in the dispersion of the shear estimates 
is a result of assuming a prior knowledge of the intrinsic elliptic¬ 
ity distribution when estimating the IA. Table [2] shows the mean 
recovered shear and IA estimates and the standard deviation of the 
estimates. 

Upon assuming a linear approximation of the F\ (|t IA |) 
function using equation (|46j, we can write a linear approximation 


of the hybrid shear estimator as 






as 


obs 

e l,i ~ 

cos (2aJ nt ) 

up?* 

obs 

e 2 ,i ~ 

sin (2aj nt )' 

up?* 


(51) 


From here we can recover an approximation of the dispersion in the 
shear estimates: 


1 

2 , 
+ 

l 

2 

/ obs 

( e 1 cos 

Vn 

2 u 2 p? t2 

u 

i 

2 , 
+ 

1 

2 

/ obs • 

(e 2 sm 

Vn 

2 u 2 p? t2 

u 




(52) 


which depends upon the correlations between the true total elliptic¬ 
ities and the true intrinsic trigonometric functions. It can be shown 
that these correlation terms can be written as 

obs / n int \ \ / obs • (n int \ \ 

e 1 cos [2a J)=(e 2 sin 12a )) 

~u + O (| 7 | | 7 IA |) + O ^ 7 IA | ) > 

(53) 


where u is a zeroth order term, which is independent of the input 
shear and IA signals, but is dependent on the form of the intrinsic 
ellipticity distribution, / (|e ran |). Hence, to first order in the shear 
and IA, the correlation terms are constant, and therefore, we can 
approximate the dispersion in the shear estimates to be 


1 

Vn 


2 i 
+ 


2u 2 P? t 


2u' 

u 


(54) 


For a Rayleigh distribution it is possible to recover the coefficient 
v! analytically, if we adopt the same assumptions used to derive 
equation (|47[>. This is found to be 


(55) 
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Therefore, we can conclude that, in the absence of shear and IA 
signals, and assuming zero measurement errors on the estimates 
of a lnt , there is a dispersion in the shear estimates which arises 
from random shape noise. With these assumptions, we found in the 
previous subsection that the dispersion in the FAO shear estimator 
was zero. For the case of the FAO estimator, a knowledge of the 
intrinsic ellipticity distribution is assumed for the random compo¬ 
nent of both the observed and intrinsic ellipticities. This allows us 
to recover estimates of the vectors (7 + 7 IA ) and 7 using only 
measurements of a obs and a mt . Using only measurements of the 
position angles eliminates the contribution of random shape noise 
in the FAO shear estimates. Flowever, the hybrid estimator requires 
measurements of e° , which contributes random shape noise to the 
estimates of the shear. This noise is, to first order, independent of 
both the shear and IA signals. 

Using equation a with the input values used in Figure[5] we 
recovered approximations for the dispersion in the shear estimates 
using the hybrid method. These were found to be = cr^ 2 = 
1.40 x 10 -2 , which are in good agreement with the values quoted 
in Table [2] 

Figure[7]shows the residual bias in the hybrid estimator. Flere 
we see a residual bias which is a result of the nonlinear propagation 
of noise on the mean trigonometric functions into estimates of the 
IA signal and hence into the shear estimates. However, we find that 
this bias is much smaller than the dispersion in the estimates in all 
of the tests we have conducted, and is negligible when we consider 
the power spectra reconstructions in Section]?] 


4 TESTS ON SIMULATIONS 

In this section we test the three estimators, described in the previous 
sections, by reconstructing the lensing and IA auto and cross-power 
spectra following the approach described in BB11. All of the simu¬ 
lated fields are assumed to be pure Gaussian fields and, as our aim 
is to demonstrate the power of the estimators to separate the shear 
and IA signals given an unbiased estimate of the intrinsic position 
angle, we ignore the effects of observational systematics. 

In all simulations we assume a ACDM background cosmol¬ 
ogy, with the matter density parameter f2 m = 0.262, the amplitude 
of density fluctuations as = 0.798, the Hubble constant Hq = 
71.4kms _1 Mpc -1 , the baryon density parameter Hb = 0.0443 
and the scalar spectral index n s = 0.962. 

We simulate the weak lensing and IA fields in three differ¬ 
ent redshift bins and include all of the possible cross-correlations 
between the fields in the different bins. The selected bins are 
0.00 < z\ < 1.40, 1.40 < 22 < 2.60 and 23 > 2.60, with the 
bin limits selected such that each bin contains approximately the 
same number density of sources. We simulate the IA signal using 
the modified non-linear alignment model introduced by |Bridle &| 
|King||2007| and we use a normalization for the IA power spectrum 
which is five times the observed SuperCOSMOS level in order to 
make it easier to see the effects we are dealing with, and we as¬ 
sume a correlation coefficient of p c = — 0.2. A detailed discussion 
of the simulated fields is given in BB 11. Here we focus on the per¬ 
formance of the estimators. 

In order to demonstrate the methods discussed, we assume that 
all of the galaxies have sufficient information to measure the IPA. 
We assume a measurement error on the a lnt (IPA) estimates with 
r.m.s. 10° and a negligible error on the ellipticity and a obs mea¬ 


surements in order to make a fair comparison between the various 
methods. 

To estimate the power spectra from the simulated observa¬ 
tions, we pixelize the sky into 3.4 x 3.4 arcmin 2 cells and assume a 
background galaxy number density of 4 arcmin -2 for each redshift 
bin; this number density is chosen to avoid the issue discussed at the 
end of Subsection |2.3| though we note that this may be achievable 
in future deep surveys with the SKA. We then reconstruct shear 
and IA maps using each of the estimators discussed in the previous 
sections. 

We estimate the recovered power spectra from the recon¬ 
structed shear and IA maps using the standard pseudo-© approach 
( |Brown & Battye|201 lb||Hivon et al.|2002| [Brown et al,|2005| . 

In the presence of noise in the shear estimates, we can write 
the general expectation value of the estimated pseudo-© power 
spectra, (5/° , as 

( cr) = C?- Y - + Cf *^ , (56) 


where the postscripts X and Y denote the fields being correlated 
and where the subscripts s and n respectively denote the signal 
and noise in that field. One can correct for biases due to noise and 
correlations between the signal and noise, such that an unbiased 
estimate of the power spectra can be recovered using a suite of 
Monte-Carlo simulations: 

cr.cf r -(cf* Y ») , 


(57) 

where the angle brackets indicate the mean over the suite of Monte- 
Carlo simulations. This is the form of the power spectra estimator 
used for the remainder of this paper. In the presence of model de¬ 
pendent noise and correlations between the signal and noise, un¬ 
biased estimates of the power spectra are only achievable if the 
Monte-Carlo simulations include the input power spectra. In a real 
analysis, this will obviously not be possible. In order to address 
this issue, we adopt an iterative approach to estimating the spectra. 
To begin with, we construct a suite of 200 Monte-Carlo simula¬ 
tions under the assumption that the input shear and IA signals are 
zero. This provides us with an initial estimate of the power spec¬ 
tra using equation As we shall see this is sufficient when us¬ 
ing the CBB and hybrid methods to recover the shear and IA esti¬ 
mates. However, it is insufficient when using the FAO estimator. It 
does, however, provide us with initial estimates of the power spec¬ 
tra. These initial estimates can then be used to construct a suite of 
improved Monte-Carlo simulations, which can be used to update 
our estimates of the power spectra. 

Figure[8]shows the reconstructed shear and IA auto and cross¬ 
power spectra for each of the three overlapping redshift bins, re¬ 
covered using the CBB estimator (black points) and using a suite 
of 200 Monte-Carlo simulations under the assumption that the in¬ 
put shear and IA signal are zero. The blue points show the recon¬ 
structed power spectra using the original BB estimator to estimate 
the shear and IA. The red curves show the input power spectra. 
From this we clearly see the success of the correction. 

Figure [9] shows the reconstructed power spectra when using 
the FAO estimator (black points). The linear form of the FAO es¬ 
timator, given in equation ([48j, has been used to reduce compu¬ 
tation time. From this we see that there is a residual bias in the 
shear power spectra which propagates into estimates of the shear- 
IA cross-power spectra. This bias is due to a dependence of the er¬ 
rors on the shear estimates on the input shear signal, as described at 
the end of Subsec tion|3. 1| This bias is not successfully corrected for 
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Figure 7. Same as for Figure[I]but for the hybrid shear estimator. Here we see a small residual bias due to the finite number of source 
galaxies. However, this bias is much smaller than the bias in the original BB estimator (Figure[lJ. 


when using noise-only Monte-Carlo simulations. However, if we 
use these estimated power spectra as the input power spectra for a 
further set of Monte-Carlo simulations, we can construct noise and 
noise-signal power spectra which include an estimate of the shear 
and IA signal. These updated noise and noise-signal power spectra 
can then be used to recover an improved estimate of the input power 
spectra. This step can be iterated until subsequent estimates of the 
power spectra are deemed consistent. The blue points in Figure [8] 
are the result of this procedure, using just one iteration. We see that 
the iterative step has, indeed, improved our estimates of the power 
spectra. 

The green points in Figure [9] show the reconstructed power 
spectra using the hybrid estimator. The linear form of the hybrid 
estimator, given in equation < |5 1 1 ), has been used to reduce com¬ 
putation time. This reconstruction did not require the use of the 
iterative procedure. It was shown, in the discussion which follows 
from equation s that the dispersion in the shear estimates is in¬ 
dependent of the input shear and IA signal to first order. Hence, the 
zero-signal noise power spectra successfully removes the noise bias 
from the power spectra without the need of the iterative procedure 
described above. 

In Figurc[l0]we show the fractional errors on the reconstructed 
power spectra. From this we see, as expected, that the errors are 
largest for the CBB estimator, and smallest for the FAO estimator. 
However, we emphasize that we have assumed a perfect knowl¬ 
edge of / (|e ran |) when using the FAO and hybrid estimators. This 
would obviously not be the case in a real analysis where uncertain¬ 
ties on our knowledge of / (|e ran |) would lead to an increase in the 
errors of the FAO and hybrid approaches. 

Table [3] shows the mean fractional bias in the shear auto and 
cross-power spectra and the IA auto and cross-power spectra recon¬ 
structions. From this we see that there is approximately an order of 
magnitude reduction in the fractional bias of the CBB estimator, as 
compared with the original BB estimator. The reduction in bias is 
generally greater when using the FAO and hybrid methods. How¬ 
ever, these methods require an accurate knowledge of /(|e ran |) 


and, for the case of the FAO estimator, an iterative method to re¬ 
move noise bias. 


5 CONCLUSIONS 

When we include a correction term into the formalism of the esti¬ 
mator introduced by BB11, we have demonstrated that the residual 
bias in the estimator, which emerges in the presence of measure¬ 
ment errors on the intrinsic position angle estimates and a non-zero 
IA signal, can be reliably reduced to negligible levels as compared 
with the original BB estimator, provided that a sufficient number 
of resolved background galaxies have reliable polarization infor¬ 
mation. When including the correction term, chance alignments of 
the measured IPA may result in substantial outliers in the distribu¬ 
tion of the shear estimates if the number of background galaxies 
is small. However, we have introduced a method which may be 
used to place constraints on the number of background galaxies 
required to recover reliable shear estimates when using this estima¬ 
tor. This restriction may require large cells, such that the number 
of source galaxies within each cell is greater than or equal to the 
minimum number of galaxies required, and hence small scale in¬ 
formation may not be attainable. 

Building upon the angle-only estimator introduced by |Whit-] 
jtaker et aL] l |2014) , we have constructed an angle-only IA estimator 
which uses IPA measurements and requires a knowledge of the in¬ 
trinsic ellipticity distribution. From here we can formulate two dis¬ 
tinct shear estimators. The first is the FAO shear estimator, which 
requires measurements of the observed position angles. The second 
is the hybrid method, which combines the angle-only IA estima¬ 
tor with the standard shear estimator and requires measurements of 
the observed ellipticities. We have demonstrated that both of these 
methods may be used to recover shear estimates which exhibit neg¬ 
ligible residual biasing as compared with the original BB estimator. 
The FAO method, however, requires the implementation of an iter¬ 
ative procedure to mitigate the effects of a signal dependent noise 
bias in the shear and shear-IA power spectra. We further emphasize 
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Figure 8. Reconstructions of the tensing and IA auto and cross-power spectra. In each panel the red curve shows the model power 
spectra. The black points show the reconstructed power spectra using the CBB estimator to estimate the shear and IA signals. The blue 
points show the reconstructions using the original BB estimator, as a comparison. From these reconstructions we clearly see that the 
residual bias has been reduced when using the CBB estimator. 
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Figure 9. Reconstructions of the lensing and IA auto and cross-power spectra. In each panel the red curve shows the model power 
spectra. The black points show the reconstructed power spectra using the FAO estimator to estimate the shear and IA signals and with 
the noise maps created under the assumption that the input shear and IA signals are zero. The blue points show the reconstructed power 
spectra recovered upon using the iterative procedure described in the main text. From this we see the success of the iterative procedure. 
The green points show the reconstructed power spectra using the hybrid estimator to estimate the shear and IA signals, with no iterative 
procedure required. 


© 2015 RAS, MNRAS 0Q0.|T|fT7l 

































































































































Separating weak leasing and intrinsic alignments 





Figure 10. The fractional error in the power spectra reconstructions shown in Figures [8] and [9] The curves show the CBB estimator 
(black), the FAO estimator (blue) and the hybrid estimator (red). Note, the input 11-13 cross-power spectrum is zero for all multipoles 
and therefore we show the error as opposed to the fractional error for that panel. 
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Spectrum 

Original BB 


Corrected BB 

FAO 

Hybrid 

Gl-Gl 

(-7.48 ±0.02) x 

10 -2 

(7.45 ± 0.29) x 10“ 3 

(1.45 ± 0.13) x 10 -3 

(0.17 ± 0.21) x 10" 3 

G2-G2 

(-3.99 ±0.10) x 

10 -3 

(9.60 ± 1.14) x 10“ 4 

(3.90 ±0.96) x 10" 4 

(1.77 ± 1.05) x 10" 4 

G3-G3 

(-1.97 ± 0.10) x 

10~ 3 

(4.02 ± 1.01) x 10“ 4 

(-6.27 ± 0.93) x 10" 4 

(0.69 ±0.98) x 10" 4 

G1-G2 

(-3.07 ± 0.01) x 

10~ 2 

(2.44 ±0.16) x 10~ 3 

(-0.38 ± 0.11) x 10" 3 

(0.26 ±0.13) x 10" 3 

G1-G3 

(-3.48 ± 0.01) x 

10 -2 

(2.68 ±0.17) x 10“ 3 

(-0.83 ±0.12) x 10" 3 

(-0.05 ±0.15) x 10" 3 

G2-G3 

(-4.73 ± 0.10) x 

10 -3 

(6.42 ± 1.03) x 10“ 4 

(-11.72 ±0.95) x 10~ 4 

(0.79 ±0.99) x 10“ 4 

11-11 

(-3.85 ± 0.01) x 

10" 1 

(3.15 ± 0.17) x 10~ 2 

(0.17 ± 0.14) x 10" 2 

(-0.27 ± 0.14) x 10" 2 

12-12 

(-3.85 ± 0.02) x 

10” 1 

(3.80 ±0.38) x 10~ 2 

(-0.12 ±0.29) x 10 -2 

(-0.39 ±0.30) x 10" 2 

13-13 

(-3.83 ± 0.09) x 

lO" 1 

(-0.66 ± 1.57) x 10" 2 

(2.79 ± 1.33) x 10" 2 

(2.31 ± 1.42) x 10“ 2 

11-12 

(-4.04 ± 0.19) x 

10" 1 

(1.15 ± 0.32) x lO -1 

(-0.56 ±0.27) x lO’ 1 

(0.33 ±0.27) x lO^ 1 

12-13 

(-4.60 ± 0.30) x 

10" 1 

(1.85 ± 0.49) x lO' 1 

(-0.09 ±0.42) x lO” 1 

(-0.23 ±0.42) x 10" 1 


Table 3. The mean fractional bias in the power spectra reconstructions across all multipoles. 


that the results presented in this paper are based on the assump¬ 
tion that the distribution, / (|e ran |), is known exactly. An incorrect 
knowledge of this distribution propagates as a multiplicative bias 
into the shear and IA estimates. However, it is expected that this 
distribution may be accurately measured using deep calibration ob¬ 
servations in future surveys. Constraints on the accuracies required, 
and the number of galaxies required to achieve these accuracies, are 
discussed in |Whittaker et al.| ( [2014j ). 

Present radio surveys, such as SuperCLASS, currently under 
observation using the JVLA and e-MERLIhQarrays, will hopefully 
provide information about the fraction of galaxies with reliable po¬ 
larization information and the expected error on the intrinsic orien¬ 
tation estimates, a a mt, provided by measurements of the PPA. This 
information is essential if we are to gain an understanding of the 
cosmological scales which may be probed using these techniques. 
In addition to this, we hope to improve our understanding of the 
impact of Faraday rotation on measurements of the PPA. It is ex¬ 
pected that this effect may be corrected for using information from 
multiple frequencies to extract the rotation measures of the source 
galaxies. 

We aim to apply these techniques to future radio surveys, such 
as with the SKA, where the number density of galaxies will be 
higher than for current radio surveys, enabling us to probe smaller 
cosmological scales. The high redshifts achieved by the SKA will 
also enable radio weak lensing to probe regions of the Universe 
which are inaccessible to other weak lensing surveys ([Brown et al. 
|2015[ >. This high redshift information will provide powerful con¬ 
straints on the evolution of large-scale structure in the Universe. 
Another exciting prospect for future radio weak lensing is the cross¬ 
correlation of radio and optical weak lensing surveys, such as the 
correlation of shear estimates from EUCLID with those from the 
SKA. This method has the advantage that the systematics in the 
two telescopes are expected to be completely uncorrelated, allow¬ 
ing the effects of systematics to be removed from shear analyses, 
while avoiding the residual effects of an incorrect calibration. 


We begin by redefining the matrix A as 

Dint N 


A = 


2/3r 

N 


E 

i=1 


WifufiJ , 


(A-l) 


and the vector b as 


b = 


2P?' 

N 


^2 Wi (e° ba ■ fi2j fn, (A-2) 


where the definition of the vector is given in equation 0 and 
where Wi is a normalized arbitrary weight assigned to each galaxy. 
In the limit N —> oo, the matrix A can be written as 


A = pt 


1 — (cos (4d mt )) — (sin (4d lnt )) 

— (sin (4d mt )) 1 + (cos (4a lnt )} 


(A-3) 


and the vector b can be written as 


b = p? 


( e obs 


>-(ef s sin (4a int )) 


} - (e° bs sin (4a int )) 


1 — cos (4d lnt ) 

(e! bB fl + cos (4d int ) 

(A-4) 

In the presence of noise on the estimates of a lnt , the trigono¬ 
metric functions above will be biased. This bias can be corrected 
for by dividing the functions by the correction term /34, as defined 
in equation ( | 14| ), such that the corrected matrix A can be written as 

pt - (cos (4d int )) - (sin (4a int )) 

1 — (sin (4d int )> /3i nt + (cos(4a mt )) 


(A-5) 


and the corrected vector b becomes 


obs 


- ( 

" (3 b 


Pt - cos (4a int ) 
pt + cos (4d int ) 


)-(e£ bs sin (4d int )} 
> - (e° bs sin (4d int )) 


(A-6) 

The CBB estimator can then be written more concisely by 
defining the matrix M, ; , where 


A/r _ ( Ft - cos (4dj nt ) 
Mi "' -sin(4d“) 


— sin (4dq nt ) 
Ft + cos (4d( nt ) 


(A-7) 


such that the final form of the estimator can be written as described 
by equations m-m- 


APPENDIX A: CORRECTING THE BROWN & BATTYE 
ESTIMATOR 

In this section we discuss the correction to the BB estimator. 

1 http://www.e-merlin.ac.uk/legacy/projects/superclass.html 
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